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ABSTRACT 

We clarify the nature of the small-diameter supernova remnant (SNR) Kes73 
and its central compact source, IE 1841—045, using X-ray data acquired with the 
ASCA Observatory. We introduce a spatio-spectral decomposition technique necessary 
to disentangle the ASCA spectrum of the compact source from the barely resolved 
shell-type remnant. The source spectrum (1-8 keV) is characterized by an absorbed 
power-law with a photon index a ~ 3.4 and ~ 3.0 x 10^^ cm“^, possibly 
non-thermal in nature. This bright X-ray source is likely a slowly spinning pulsar, 
whose detection is reported in our companion paper (Vasisht & Gotthelf 1997). The 
SNR spectrum is characteristic of a thermal plasma, with kT ~ 0.7 keV, and emission 
lines typical of a young remnants. The element Mg and possibly O and Ne are found 
to be over-abundant, qualitatively suggesting an origin from a massive progenitor. 

We hnd that Kes73 is a young ( 2000 yr) type Il/Ib SNR containing a neutron 

star pulsar spinning anomalously slow for its age. Kes73 is yet another member of a 
growing class of SNRs containing radio-quiet compact sources with a hard spectral 
signature. 


Subject headings: stars: individual (Kes73, IE 1841—045) — stars: neutron — 
supernova remnants — X-rays: stars 


1. Introduction 

The supernova remnant Kes73 is a limb - brightened, shell-type radio remnant ~ 4' in 
diameter, located along the Galactic plane (G 27.4-1-0.0). The X-ray remnant is comparable in 
size to the bright radio shell, but is dominated by localized emission knots, interpreted as clumpy, 
diffuse emission in the interior, possibly due to fluorescence from reverse shock (Helfand et al. 
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1994 using ROSAT HRI images). Centered on the remnant is the unresolved compact X-ray 
source, IE 1841—045, discovered with the Einstein Observatory (Kriss et ah 1985). Like the 
supernova remnants RCW 103 and CTB 109, the radio observations show no evidence for a 
plerionic component or a localized radio counterpart to IE 1841—045, with a flux limit Egcm < 0.6 
mJy (Kriss et al. 1985). An HI absorption-based distance determination shows that Kes73 lies 
~ 7.0 kpc away (Sanbonmatsu & Helfand 1992). 

The nature of Kes 73 and in particular, its compact source, has been considered extensively 
in two earlier studies. Kriss et al. (1985) were unable to distinguish between thermal and 
non-thermal emission from IE 1841—045; they suggested as the origin of the emission either a hot 
thermal neutron star (NS), a crab-like pulsar, a plerionic nebula, or an accretion powered binary. 
A more recent study by Helfand et al. (1994) reconsidered these possibilities using deep ROSAT 
HRI pointings. They inferred a hard spectral component for IE 1841—045, placed a limit on its 
long term variability of a factor of two change in flux, and suggested an accreting binary origin for 
the central source. 

In this Letter, we use the broad-band ASCA X-ray observations of Kes 73 to further clarify 
its nature. We present a novel technique for isolating the background subtracted spectrum of 
the central source. We show that the spectrum of IE 1841—045 is represented by a steep power 
law, possibly of non-thermal origin, and suggest that Kes 73 is a young ( 2000 yrs) Type H/Ib 

SNR, the product of a massive progenitor. In our companion paper (Vasisht & Gotthelf, referred 
herein as GV97) we report our discovery of a 11.8 s periodicity from IE 1841—045, likely from an 
anomalous pulsar, the stellar remnant of the supernova explosion. The date and typing of Kes 73 
are critical to understand the nature of the pulsar. We present Kes 73 as another member of an 
emerging class of young thermal remnants which are found to contain a radio-quiet, hard X-ray 
point source. 


2. Observations 

ASCA (Tanaka et al. 1994) conducted a two day observation of Kes 73 on 11-12 Oct 1993. 
In this study, we use data acquired with the two Solid State Imaging Spectrometers (SIS-0 and 
SIS-1), made available from the public archive. These sensors are sensitive to X-rays in the 
0.4 — 10.0 keV band, with a nominal spectral resolution of 2% at 6 keV ( E ^/^). The spatial 
resolution is limited by the X-ray mirrors, whose azimuthally averaged point spread function (PSE) 
is characterized by a narrow core of EWHM 50'', and extended wings that result in a half-power 
diameter of 3' (Jalota et al. 1993). The SIS data were acquired in 1-CCD mode, read out every 
4 sec, using a combination of EAINT and BRIGHT data modes. The target was centered on the 
IT X 11' field-of-view of SIS-0 CCD-1. All data were edited using the standard REVl screening 
criteria which resulted in an effective exposure of ~ 34 ksecs per sensor. A description of the data 
acquired with the other two focal plane instruments is presented in our companion paper (GV97). 
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The broad-band image of Kes 73 from the combined SIS cameras contains 93 kcounts, with 
a count rate of 1.38 cps. In Plate 1 we display images of the Kes 73 region in the soft- and 
hard-band, below and above 2.5 keV, respectively. The images are centered on the point-like 
source, whose position is consistent with that of IE 1841—045, to within the formal SIS error circle 
(Gotthelf 1996). The hard-band image clearly shows the clover-leaf shaped PSF, consistent with 
an unresolved point source, whereas the soft-band image suggest an additional diffuse structure. 
The FWHM of the normalized brightness distribution profiles in the two bands are 44'' (hard) and 
108" (soft) (see Fig 1). This difference is more apparent in the deconvolved images in the lower 
panel of Plate 1. These clearly show that the thermal X-ray shell is contained in the soft-band and 
surrounds the harder, bright compact source. The images have been restored with 50 iterations of 
the Lucy algorithm (Lucy 1974) to deconvolve the complex shape of the PSF, and are displayed 
in an identical fashion. 

An image of the SNR nebulosity is made by subtracting the compact component from the 
soft-band image. For this, we used as a model PSF the SIS image of the moderately bright 
point-like source EX Hydra, normalized by the derived source spectrum (see below) to estimate the 
relative contribution of the compact source to the total emission below 2.5 keV. The subtracted 
SIS image closely follows the brightness distribution in the smoothed Rosat HRI image (Plate Ic), 
reproducing the three areas of enhanced emission. Other features are easily ascribed to differences 
in the spectral response of the two instruments. See Hwang & Gotthelf (1997) for a discussion on 
interpreting features in the processed images. 


3. Isolating the compact source spectrum 

The broad winged ASGA PSF does not allow spatial isolation of the point-source emission 
from that of the SNR. Instead we adopt a spatio-spectral decomposition method detailed in Wang 
& Gotthelf (1997) (referred herein as WG97) to separate the spectra of the two components. 

We then use the decomposed point-source spectrum to: (i) verify and consolidate the fit to the 
point-source spectrum in the total spectrum of the Kes 73 region (see §4.2) and (ii) compute fluxes 
and luminosities of the point-source. Below we outline a modified version of the WG97 method, 
germane to a point source embedded within a SNR shell. 

We start with the assumption that the spatial distribution of the Kes 73 emission consists 
of two components, i.e., a point-source embedded in a diffuse nebula. We then decompose the 
source spectrum from the underlying nebula (-|- background), by simultaneously solving for the 
source and nebular spectra using the ratio of observed and expected counts in concentric regions; 
we search for deviations from the radial-average profile centered on the source, from that expected 
for a point-source. 

The original method, discussed in WG97, assumes a uniform background. From the Rosat 
HRI morphology, we know this is not the case for Kes 73, for which the background consists of both 
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the thermal-shell emission and the field-background. Therefore, we use the HRI data to estimate 
a correction factor to a uniform background from the relative counts in the source and nebular 
region. Our procedure is as follows; (i) We first extract spectra from two annuli centered on the 
source (see Fig. 1). The first annulus is a small circle (r = \'.2) encompassing the central pulsar 
counts. The second annulus (lf5 > r > 2'.7) is chosen to enclose the bulk of the shell emission, (ii) 
Next, we compute the expected radial profile for a point-source in these regions from a similar 
1-CCD mode observation of EX Hydra. To compensate for the non-uniform SNR distribution, 
we compute a nebular correction factor using the relative counts per pixel between the center 
and edge of the projected HRI X-ray nebula (hi) With the above information, we decompose the 
spectra in the two annuli into a source spectrum and a nebular spectrum using eqns. A2 3-4 of 
WG97, applied to each spectral channel. 

The spectra separated into two unmixed components: a line dominated spectra expected 
from thermal emission of a shocked SNR shell and a steep power-law continuum spectrum for the 
compact source. These are shown on Plate la &: lb. No coercion or prejudice is used to force the 
clean separation into the distinct nebula and source spectra. The background subtracted source 
spectrum (Plate la) accounts for the harder (> 2.5 keV) emission and provides an independent 
absorption measurement. The spectrum is trivially fitted with an absorbed power-law with a 
photon index T = 3.4 ± 0.3 and Nh — 3.0 ± 0.4 x 10^^ cm“^. The inferred unabsorbed luminosity 
is Lx(l-0 — 10.0 keV) ~ 3 x 10^® df erg s“^, for an assumed distance of 7.0 di kpc (see Table 
1). The spectral shape (oc E~‘^) is in accordance with the steep power-law spectra seen in other 
anomalous X-ray pulsars (see Corbet et al. 1995). We now use this fit to constrain the power-law 
emission component in the combined fit to the Kes73 spectrum. 


4. The SNR nebula spectrum 

The SIS energy spectrum photons were selected from within a circular emission region of 
radius ~ 4', limited by the size of the CCD. The broad ASCA PSF makes background estimation 
extremely difficult, since the source flux from a diffuse object extends over most of the CCD chip. 
For this analysis, we extract a background spectrum from nearby archival pointings of the Calactic 
ridge. The resulting background-subtracted spectrum of the shell plus compact source is shown 
in Fig 2. The line-dominated spectrum, typical of those seen from SNRs, suffers high foreground 
absorption in the energy range below ~ 1.5 keV. In the energy range spanning ~ 1 — 10 keV, we 
see K-shell emission from highly ionized ions of Mg, Si, S, Ca, and Ar. The identified emission 
lines and their characteristic parameters are displayed in Table 2. 

Collisional ionization equilibrium (Raymond & Smith 1977; Mewe et al. 1985 and references 
therein) and thermal bremsstrahlung models with Caussians for line emission, are fit to the 
spectra. There are large residuals, in either case, in the hard-band (2.5 - 10.0 keV) suggesting an 
extra emission component. The image analysis, spanning that energy range, clearly demonstrates 
that this component is mostly emission from the compact source, for which we add a single 
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power-law to the fit. A bremsstrahlung continuum, with Gaussians and a power-law tail provides 
the best fit (see Table 1). In contrast, the single temperature Raymond-Smith model with 
fixed relative metal abundances and a power-law, is a relatively poor fit; it shows large negative 
residuals mainly near the over-abundant Mg feature, and a relatively large metal abundance of 
~ 2.0 cosmic, driven mainly by the strong Si (1.83 keV) feature. In addition, we note excess fiux 
in all our fits, at the lower energy range, 0.5 — 0.9 keV, which we tentatively ascribe to O and Ne. 

5. Discussion 

A lower limit on the SNR age can be derived assuming free expansion of the spherical 
remnant. For a typical maximum velocity of 10^ km s“^ seen in Type II SN, the projected size of 
Kes73, Rs ~ 4.7 dj pc, contrains its age to be Tg < 470 yrs. Since most of the thermal emission 
is due to shocked matter of electron temperature kT^ = 0.8 keV, we compute a shock speed 
Vs = (16kTe/3^mp)^/2 900 km s This assumes ion and electron equilibrium behind the 

shock {p, = 0.6 for cosmic abundances). Inefficient electron heating would result in a larger shock 
velocity, i.e., Vg ^ 900 km s“^. If we assume that the remnant has entered a well developed Sedov 
phase, Rg = 2.5 vsTs, and therefore Tg ^ 2.2 x 10^ yr. This age is further reduced if the remnant is 
not fully in Sedov phase. From spectral fitting, the total thermal (1 — 10 keV) luminosity of the 
SNR is Lx ~ 3 X erg s“^ (this estimate excludes the uncertain, highly-absorbed emission 

from O and Ne). The instantaneous power radiated from a shell is Lg{t) = {16iT/3)Rl{t)n‘^A{T), 
where Uq is the mean pre-shock particle density and A(T) = 1.0 x -|- 2.3 x 10“^^Tg'® 

erg cm^ s“^ is the cooling function (McCray 1987), and Tq is the kTe in units of 10® K. We derive 
an Uo ~ 0.8^7 cm“^. Under the strong-shock assumption the post-shock electron density is ~ 3 

cm“^. 

Kes 73 is evidently young, with the shocked plasma still ionizing and thus non-equilibrium 
effects in the ionization balance are important. We determined the diagnostic parameters for the 
non-equilibrium ionization (NEI) plasma from the line intensity ratios of He-like-Ka to H-like-Ka, 
as well as, He-like-K/3 to He-like-Ka ionic transitions of individual ions. Si and S in this case. 
The former ratio is a measure of the ionization degree and is dependent upon both the electron 
temperature kTe and the ionization age n-et, the product of the electron density and time since 
the gas was last heated by the shock (Itoh 1977). The latter ratio, however, is solely a function of 
kTf,. The obtained He-hke-K/3 to He-like-Ka ratio for Si and S is 0.087 ± 0.015 and 0.089T0.039, 
respectively. The He-like-Ka to H-like-Ka intensity ratio of Si and S are 22 ± 9 and 15 ± 9, 
respectively. We measure a kTe ~ 0.75 — 0.90 keV and 0.71 — 0.97 keV for Si and S, respectively, 
assuming a single uniform plasma continuum and a power-law. These are in accordance with the 
kTe derived from the single temperature bremsstrahlung continuum. The ionization parameter was 
estimated to lie in the range Uet ~ 11.0 — 11.4 (Masai 1984). Using the post-shock density inferred 
via normalization to the continuum component, we get t ~ 1800 yr. Several independent 
arguments, therefore, suggest that Kes 73 is a young SNR, ~ 2000 yr-old, as indicated by earlier 



- 6 - 


studies (Helfand et al. 1994). 

We estimate the total mass of the swept-up interstellar gas to be Mg ~ 8.8dj M©. This 
corresponds to typical envelop masses ejected in a Type II supernova (progenitor mass M > 8 
Mq), suggesting that the SNR dynamics could well still be in a transitional stage between the free 
expansion and Sedov phases. This notion is qualitatively consistent with the development of a 
strong reverse shock that can heat the metal-rich gas and cause the diffuse emission to be seen in 
the SNR interior. Additionally, the detection of strong Mg, Si, S, and Ar emission and perhaps 
accompanying O and Ne emission suggests ejecta-dominated gas. During their evolution, massive 
stars are expected to produce large amounts of 0-group elements which are ejected during the 
supernova (Thieleman, Nomoto & Hashimoto 1994; Woosley 1991). These patterns have been 
observed in the X-ray spectra of 0-rich SNR (see Hughes & Singh 1994; Hayashi et al. 1994). 
Here the emission is highly absorbed and does not permit quantitative estimates of O and Ne 
emission (although, upper-limits and observed photon excess in the spectral range 0.4 - 1 keV are 
consistent with abundant O and Ne). 

The thermal spectrum of Kes 73 is remarkably similar to that of two young, distant (therefore, 
high-absorbed) Type H/Ib SNRs, G 11.2—0.3 and ROW 103. The former is a historical SNR 
(supernova A.D. 386) and hence approximately the same age as Kes 73 (Vasisht et al. 1996). 
Unlike Kes 73, it contains a weak, but extended hard X-ray plerionic core. ROW 103 is a virtual 
twin of Kes 73; both have radio-quiet, point-like X-ray sources in their centers, share morphological 
similarities, while neither shows evidence for a radio or X-ray plerion (see Gotthelf et al. 1997 
and refs therein). We infer a minimum energy in particles and nebular magnetic fields in the Kes 
73 core to be Emin < 10^^ erg, assuming Crab-like parameters and equipartition between 

magnetic fields and relativistic particles (Pacholczyk 1970). 

The lack of observed plerionic emission can be qualitatively explained in the following manner. 
If the Kes73 pulsar has a large dipolar field B ~ 10^^ G (VG97), a weak plerion is a natural 
consequence at an age of ~ 2 x 10^ yr. Such a pulsar loses most of its initial spin energy in 
a matter of ~ years. The released pulsar wind would immediately suffer strong adiabatic and 
synchrotron losses. This leads to a bright plerion at an early age (~ 10 — 100 yrs) with subsequent 
rapid decline in surface brightness. This effect is shown by Bhattacharya (1990) in his paper on 
the morphology of SNR with central pulsars. He shows that for pulsars with R-fields spanning the 
range 10^^ G to 1.5 x 10^^ G, the ones with the highest B have the faintest plerions at an elapsed 
time t ~ 10^ yr. For magnetars, B ~ 10^^ G, this rapid decline in surface brightness may be much 
more pronounced. 
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Fig. 1. " Radial profiles of Kes 73 in the hard- and soft-bands (above and below 2.5 keV). 

The hard-band profile is consistent with an unresolved point source, while the soft-band profile 
indicates additional diffuse emission. The dotted vertical lines indicate the regions used in the 
spectral decomposition (see §3). 

Fig. 2. ~ The SIS spectrum and model of the Kes 73 region: The upper panel represents a 
background subtracted pulse-height spectrum of the supernova remnant in the entire region 
(r ^ 4'), and the best model function of a thermal plasma consisting of a bremsstrahlung 
continuum including Gaussians emission lines (see table 2), and a power-law component to 
represent the compact source emission. SIS-0 and SIS-1 data have been summed. Crosses denote 
the actual data and the histogram indicates the model. Residuals are shown in the lower panel 
with la error bars (See tables 1 & 2 for the fitted parameters). 

Plate 1 - Images of Kes 73 using data from both SIS cameras. The brightness distribution is 
exposure corrected and smoothed. The plots are centered on the peak (broad-band) emission and 
displayed on a linear scale with contours in 16 uniform increments, scaled to the image max. Top 
Left - SIS image of Kes 73 in the hard-band (2.5 - 10.0 keV); the image is consistend with an 
unresolved SIS point source. Top Right - SIS image of Kes 73 in the soft-band (0.5 - 2.5 keV); 
evident is a additional strong diffuse component. Bottom Left - The above hard-band image 
deconvolved using the Lucy restoration method. Bottom Right - The above soft-band image 
deconvolved using the Lucy restoration method. 

Plate 2 - Decomposed SIS spectra of Kes 73. (Upper panels) The computed background 
subtracted spectrum (crosses) of the central compact source with a single absorbed power law fit 
(histogram). The best fit value for the photon index is 3.4. (Bottom panels) The computed 
spectrum for the thermal-nebula -|- background (crosses), fit with a kT = 0.6 keV Raymond-Smith 
thermal plasma model (histogram); residuals are apparent around O, Ne, and Mg. In the energy 
band below about 2.5 keV, the thermal component dominates, whereas the central power-law 
component is dominant in the harder band. The plots are shown scaled identically for ease of 
comparison. 

Plate 2 (bottom). The Rosat HRI brightness image is overlayed with the SIS contours of the 
thermal emission from Kes 73 after subtracting off the non-thermal component. The method used 
to create this image is described in the text. The HRI image is smoothed with a boxcar filter to 
bring out features correlated on the same spatial scales as the contoured image. 



Table 1. Spectral Models and Fit Parameters 


Model 

Continuum 
(F; kT) 

Nh 

(10^^ cm“^) 

xl 

— The 

Power-law 

Compact Source Spectrum — 

3.1-3.7 2.7-3.4 

1.0(151) 

Bremsstralung 

1.7-2.2 

1.9-2.4 

1.0(151) 

Blackbody 

0.6-0.7 

1.0- 1.4 

1.0(149) 

+ power-law 

0.9-3.6 

- 

- 


— The Source + Nebular Spectrum — 


Bremss -I- lines 

0.5-0.9 

2.1 - 2.2 

1.0(187) 

+ power-law 

3.4 fixed 

- 

- 

R&S 

0.6 

2.9 

3.2(207) 

+ power-law 

3.4 fixed 

- 

- 

Decomposed Nebular Spectrum — 


Bremss -I- lines 

0.6-0.8 

1.6-2.3 

1.0(74) 

R&S 

0.57-0.62 

2.7- 2.9 

3.2(88) 


Fits using ASCA SIS data between 0.9 — 8.0 keV. All 
continuum fits values are quoted in units of keV; all 
power-laws slopes are given as photon indexes iT oc 
Raymond & Smith (R&S) fits with abundance set to 2 
times cosmic. Nebular fits include 11 lines of Mg, Si, S, 
Ar, Ca. Gaussian line fits required ct = 20 eV. Errors 
are the formal 90% confidence limit for one interesting 
parameter. 



Table 2. 


Measured Emission Lines in Kes 73 


Line 

Name 

Lab 

Energy 

Line 

Center 

Norm 

10“^ cts/cm^/s 

0 Hen. 

0.65 

0.6±0.1 

< 0.2 

Ne Hea.... 

0.92 

0.9±0.1 

<5 

Fe-L blends 



< 0.1 

Mg Hea ... 

1.35 

1.340±0.003 

4.0 ±0.5 

Mg Lya ... 

1.48 

1.47±0.05 

0.14 ±0.14 

Si Hea. 

- 1.86 

1.848±0.003 

3.1 ±0.3 

Si Lya. 

2.00 

1.97±0.02 

0.14 ±0.06 

Si He3p-l-4p 

~ 2.20 

2.20L0.01 

0.27 ±0.04 

Si Ly/3. 

2.38 

2.37±0.03 

0.25 ±0.04 

S Hea. 

~ 2.45 

2.457±0.005 

0.78 ±0.06 

S Lya. 

2.62 

2.64±0.14 

0.05 ±0.03 

S Hed. 

2.89 

2.89±0.03 

0.07 ±0.03 

Ar Hea.... 

~ 3.1 

3.08±0.03 

0.11 ±0.02 

Ar/Ca Hea 

~ 3.8 

3.8±5.0 

0.01 ±0.01 


A single temperature bremsstrahlung with power-law 
model is assumed for the continuum emission. All energies 
are quoted in keV. Errors are the formal 90% confidence 
limit for one interesting parameter. Line normalizations 
are given as total photons per cm^ per sec in the line. 














Source 

circle 



1 



RADIUS (ARCMIN) 



2 

CHANNEL (KEV) 


5 


evg 29-May-1997 15:29 





















